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Abstract Using Hinode EUV Imaging Spectrometer (EIS) spectra recorded 
daily at Sun center from the end of 2006 to early 2011, we studied the long- 
term evolution of the quiet corona. The light curves of the higher temperature 
emission lines exhibit larger variations in sync with the solar activity cycle 
while the cooler lines show reduced modulation. Our study shows that the high 
temperature component of the corona changes in quiet regions, even though the 
coronal electron density remains almost constant there. The results suggest that 
heat input to the quiet corona varies with the solar activity cycle. 
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1. Introduction 



Rad iation from the c orona changes dramatically over the solar activity cy- 
cle. Peres ei"all ( 2000h stud ied the long-term va riation of the soft X-ray flux 
measured by Yohkoh/SXT ()Tsuneta et all Il99ll) and found that the emission 
measure of the whole Sun changes by orders of magnitude between the solar 
activity maximum and the minimum. The irradiances of EUV coronal emission 
lines exhibit a large variatio n in sync with the solar cycle, wi th larger amplitude 
at the higher temperatures |Del Zanna and Andrettal . 120111) . Since the coronal 
irradiance is thought to be a measure of the magnetic activity of the whole 
corona, it is important to understand what features on the Sun cause the irradi- 
ance variation. The variation of the irradiance is mainly attributed to magnetic 
field concentrations associated with active regions, fe.Q.. lFoukaill990[) . however, 
Schiihle et al. I (|2nnnh demonstrated that UV emission lines from the lower corona 
increased in quiet regions during the ascending phase of the solar activity cycle. 
Their results suggest that magnetic activity in quiet regions also changes with 
the solar cycle. 
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To quantify the variation in quiet regions, coronal emission lines formed at 
high temperatures need to be studied. We investigate the long-term variation 
of the EUV radian ce recorded daily at Sun center b y the EUV Imaging Spec- 



trometer fEIS: .Culhane et all 120071 ) aboard Hinode (jKosugi et g/J . 120071 ). The 



paper is organized as follows; the data reduction and calibration procedure are 
described in Section [5J The long-term variation of emission lines and electron 
density estimation are presented in Section [3] The results and interpretations 
are discussed in Section H) 



2. Data Reduction 

The EIS SYNOPOOl and SYNOP005 studies are regular observational programs 
recording EUV spectra near Sun center on a daily basis. Since the Sun is rotating 
at a rate of 28 days near the equator, the time series smoothed over the rotation 
period provides average properties of the corona near Sun center. However, one 
of the challenges in determining a long-term variation of the corona is that the 
sensitivity of EIS gradually decreases with time. 

A detailed study of the EIS sensitivity change was carried out by using the 
He II A256.32 A emission line recorded using the regular quiet region observation 
studies SYNOP002 and SYNOP006 near disk center (Mariska 2012, private 
communication). Those data suggest that the degradation of the EIS sensitivity 
can be expressed as 

e = {exp(i/ri)-f exp(i/r2)}/2, (1) 

where t denotes the time from Hinode launch on 22 September 2006, and the 
e-folding times ti and T2 are 467 days and 11311 days, respectively. In this study, 
the degradation of the EIS sensitivity is compensated for using Equation (1). 

A large number of individual emission line profile fits go into each averaged 
data point used to compute the sensitivity correction, and there are many possi- 
ble ways to estimate the error in each averaged data point. Two obvious choices 
are to weight each averaged data point by the number of individual measure- 
ments that went into its determination or to weight each point by the standard 
deviation of all the data points that went into the average. Each weighting results 
in nearly the same results for the initial intensity and the time constant for the 
shorter decay time, but different answers for the longer decay constant. The 
errors for the fits, however, are nearly the same: less than 5% for the initial 
intensity, about 15% for the shorter decay constant, and about 40% for the 
longer decay constant. The error for the longer decay constant is larger because 
the value itself is large — several decades. Thus the trend is very small over the 
time interval of the data fitted in this study. 

The radiance of the spectral lines listed in Table [1] is deduced by fitting a 
Gaussian function to the spectrum integrated over 100 arcsec in the north-south 
direction along the slit. This integration along the slit is intended to improve the 
signal to noise ratio of the spectrum and to smooth out small-scale structures, 
such as coronal bright points, along the slit. The spatial extent in east-west 
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Figure 1. Time series of X-ray flux along the central meridian. Dashed lines in the middle 
mark the region studied with EIS. 



direction is 1 arcsec since the spectrum is recorded with 1 arcsec wide sht. These 
emission hnes exhibit sufficient counts in quiet regions for a useful analysis. After 
correcting for the long-term degradation of the EIS sensitivity denoted by Equa- 
tion (1), data number (DN) is converted to physical units (e rg cm~^ s~^ sr" 



by using the pre-flig ht calibration of the EIS effective area (jLang et all 12006 



ICulhane et al\ . 2007 ). which is provided in the EIS branch of Solar Software 
(SSW). 

The light curve for each emission line is further smoothed along the tem- 
poral axis by convolving it with a Gaussian function with a full width at half 
maximum (FWHM) of 28 days — the solar rotation period. This is intended to 
suppress the radiance variation caused by the app earance and disappearance o f 
any small active regions during a solar rotation (IWoods and Rottmanl . 2005h . 
The smoothing process at a time to is performed by convolving with a Gaussian 
weighting function 



where FWHM = 2^2 log 2a. 



3. Results 

3.1. Location of Active Regions 

To make sure that the EIS observations at Sun center are not directly influ- 
enced by active regions, we constructed a time-space map by lining up stripes of 
central meridian d ata in daily images obtained by the X-Ray Telescope (XRT; 



Golub et all 120071) on Hinode. Figure [T] shows a time-space map composed of 
data obtained with either the ALpoly or ALmesh filters, which are available for 
most of the period. Sinc e the tempera ture responses of these two filters are almost 



identical above 2 MK ( Golub et [2007), the appearance of active regions in 



these filters should be quite similar. The time-space map is smoothed along the 
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Table 1. Fitted radiance at the beginning of Hinode ob- 
servations (December 2006) and at the minimum (February 
2009). 



Ion 


A 


logTe 


Dec 2006 


Feb 2009 


Ratio 




[A] 


[K] 


[erg cm~^ 


s-i sr-i] 




Fe XIV ^ 


264.79 


6.25 


38 


7.8 


4.8 


Fe XIII 


202.04 


6.20 


207 


33 


6.2 


Fe XII 


193.51 


6.15 


225 


70 


3.2 


Si X 


258.37 


6.15 


75 


35 


2.1 


Fe XI 


180.40 


6.10 


426 


282 


1.5 


Fe X 


184.54 


6.00 


103 


118 


0.86 


Si VII 


275.36 


5.80 


16 


23 


0.71 


Fe VIII 


185.21 


5.65 


31 


44 


0.72 


He II 


256.32 


4.85 


256 


229 


1.1 



^With a blending Fe XI 264.77A 



temporal axis by convolving it with the Gaussian weighting function denoted by 
Equation (2). The pixel size in north-south direction remains at 2 arcsec. 

The region recorded by the EIS observations is indicated with horizontal 
dashed lines. Active regions passing the meridian appear as bright features in 
the plot. Note that the EIS data are always recorded at Sun center viewed from 
the Earth, and the tilt of the solar rotation axis is not compensated for. Active 
regions occasionally appeared near Sun center from 2007 to the middle of 2008, 
while they are almost absent in 2009. In 2010, they again started to appear 
regularly, but distant from Sun center. The plot demonstrates that with only 
one exception in 2008 the EIS observing region is not covered by active regions. 
Therefore, the radiances deduced from the EIS observations are not directly 
affected by an increase or a decrease in the number of active regions. 

3.2. EUV Radiance Variation 

Figure [5] shows the smoothed light curve in each emission line. The measured 
radiance is smoothed by convolving with a Gaussian function with a width of 
the 28-day solar rotation period (see Section Since the light curves are nearly 
symmetric with respect to the minimum in 2009, the smoothed radiance is fitted 
with a second-order polynomial function to determine the long-term variation. 
The fitted polynomial functions are overplotted using dotted curves in Figure 
m and show a reasonable agreement with the light curves. In the following, 
the radiance in each emission line is determined by using the fitted polynomial 
functions. The lowest radiance for Fe xii is attained in February 2009, which is 
regarded as the minimum of the long-term variation. (The NOAA Space Weather 
Prediction Center lists the solar cycle minimum date as December 2008.) The 
radiance values at the beginning of Hinode observations (December 2006) and 
at the minimum (February 2009) are summarized in Table [TJ 

Emission lines hotter than Fe x tend to show larger variation ratios with 
increasing temperature. The only exception is the Fe xiv line, which has a smaller 
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Figure 2. Radiances measured in different emission lines are plotted as solid curves. The 
dotted curve in each panel show a second-order polynomial fit to the measured radiance. 



ratio than that of the Fe xiii hne. A possible effect of a blending emission line 
is discussed in Section [531 On the contrary, the low-temperature emission lines 
between Fe viii and Fe X show constant or a slightly larger radiance at the 
minimum. Although the sensitivity change is determined from the He ii emission 
line alone, a wavelength dependent degradation of the instrument is not likely 
the reason for this trend, because Fe viii and Si vii are observed in different EIS 
wavelength channels and yet exhibit a similar variation (Table [T]) . 



3.3. Spatial Distribution of the Radiance 



To study the spatial distribution of the radiance. Figure [3] compares the his- 
togram of the Fe Xii A193.5lA radiance in each pixel of EIS recorded in De- 
cember 2006 and in February 2009. The intense Fe Xii A193.5lA emission line 
allows us to analyze the radiance in an individual pixel (1" xl") without sum- 
ming, which can differentiate between localized and dispersed brightenings. The 
histogram is peaked at 200 erg cm~^ s~^ sr~^ in the December 2006 data (solid 
line), while that in February 2009 is at 60 erg cm^^ s^^ sr^^ (dotted line). The 
ratio of these values is close to the mean radiance ratio of 3.2 in Table [T] The 
histogram demonstrates that the long-term variation of the average radiance is 
not caused by an increase or a decrease of localized bright features, but instead 
by a variation of the entire region. 



SOLA: ms.tex; 3 March 2013; 9:10; p. 5 



S. Kamio and J.T. Mariska 




Radiance [erg cm"'' s"' sr"'] 

Figure 3. Histogram of the Fe xn A193.5lA radiance. Solid and dotted lines show the 
distributions in the December 2006 data and in the February 2009 data, respectively. 
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Figure 4. Time history of the electron density derived from the Si X line pair. 



3.4. Electron Density 

The electron density is estimated from the line ratio of Si x A258.37 A/261.04 A. 
The line ratio is deduced from the Si X radiance integrated over 100 arcsec 
along the slit. The relationship between electron density an d the line ratio is 



calculated using v ersion 7.1 of the CHIANTI atomic database (|Dere et a/.l . ll997 



Landi et aLl . l2012r ). FigureHlpresents the electron density variation. The obtained 
electron density is smoothed along the temporal axis in the same manner as 
the radiance in Section 13.21 The dotted line plotted over Figure |4] indicates 
a constant density of 3.8 x 10* cm~^. In contrast to the radiance variation 
in the Si x emission line, the electron density shows no clear trend of solar 
cycle variation. The uncertainly of the radiance due to the spectrum fitting is 
estimated to be smaller than 10%, which corresponds to an uncertainty of 30% 
in the density estimation. Note that the line ratio analysis depends solely on the 
relative radiance and is not affected by the uncertainty of the absolute radiance. 
In the case of an isothermal plasma, the observed radiance / is denoted as 

I ^ fdnlg{T,n,) (3) 
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Figure 5. Results from DEM analysis in December 2006 (left panel) and February 2009 (right 
panel) . 



where /, d, and g{T, Uc) are the fiUing factor of the emitting plasma, the depth 
over which the emission Hne is formed, and the contribution function to the Hne, 
respectively. To examine whether the density variation is responsible for the 
radiance variation, a square root of the Si x A258.37A radiance is overplotted 
in Figure |4] with a dashed curve. If the density variation is the primary cause 
of radiance variation, the density should agree with the square root of radiance. 
Figure H] indicates that the density variation is not sufhcient to reproduce the 
observed radiance variation. Other factors, the filling factor /, depth of the 
corona d, or the contribution function, are thus responsible for the radiance 
variation. 



3.5. DEM Analysis 

A differential emission measure (DEM) analysis is carried out with the aid of 
the CHIANTI atomic database. All emission lines in Table [T] except for He ii are 
used in the analysis. Since a sufficient number of emission lines are available only 
within a temperature range of log To = 5.65-6.25, we restrict the analysis to that 
interval. A DEM in that temperature range is computed by using the procedure 
provided in the CHIANTI atomic database. It deduces an optimum DEM so 
that the DEM multiplied by the contribution function of each emission line 
reproduces the measured radiance. Figure [S] presents the DEM at two different 
times. In December 2006, the DEM peak is located at logTc = 6.05, while 
that at the minimum in February 2009 the peak is found at 5.95. The high 
temperature tail of the DEM is extended higher in the former period than the 
minimum period. This explains the trend that the higher temperature emission 
lines exhibit greater decreases in the minimum period (Figure[2). In contrast, the 
lower temperature emission lines (Fe viii and Si vii) indicate slight increases in 
the minimum period as the DEM peak is shifted towards the lower temperature. 
These results show that the properties of the quiet region corona change with 
the solar activity cycle. 

The radiance in the Fe xiv A264.79 A line is a special case. It turns out that 
the radiance of Fc xiv A264.79 A is overtaken by a blending Fe xi A264.77 A 
feature during the minimum period. According to the CHIANTI estimation of 
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the radiance, Fe xiv is five times larger than Fe xi in December 2006, while 
Fe XIV gets smaller than Fe xi at the minimum in February 2009. Because of the 
Fe XI A264.77 A blending, the measured radiance at A264.79 A (Table[T]) does not 
decrease much, resulting in a smaller variation than the Fe XIII A202.04 A light 
curve. Thus care should be taken in analyzing the Fe xiv A264.79 A radiance. 



4. Discussion 

We found a long-term variation of the EUV emission lines in quiet regions in sync 
with the solar activity cycle. The time series of the EIS spectra was recorded at 
Sun center, where there is minimal contamination from active regions (Figure[T]). 
The light curves of emission lines above 1 MK clearly show that the quiet region 
corona varies with the solar activity cycle (Figure [5]) . The radiance variation is 
attributed to the coronal temperature rather than the coronal electron density. 
As shown in Figure HI the observed density variation is inadequate to reproduce 
the observed radiance variation. The temperature dependent variation of the 
light curves is explained as a modulation of the high temperature component in 
th e corona (Figure El). 



Orlando. Peres, and Reald tOOlh claimed that the average temperature of 



low X-ray flux regions observed by the Yohkoh/SXT remained constant around 
1.5 MK between the solar maximum and the minimum. A possible explanation 
for this disagreement is that the sensitivity of SXT is low below 1 MK, and hence 
is not optimum for diagnosing the corona around 1 MK. The EIS observations 
provided us with several emission lines at that temperature to allow a more 
detailed study. 

Our study proves that the corona in quiet regions is not constant, but changes 
with the solar activity cycle. It suggests that the heat input to the quiet corona 
changes with the solar activity cycle. In addition, the spatial distribution of 
the radiance indicates that the radiance variation in the quiet region is a dis- 
persed phenomenon rather than localized features (Figure [S]). Since a nearly 
linear relationship has been established between tota l magnetic flux and X-ray 



radiance in quiet regions, active regions, and stars (jPevtsov et all 120031) . our 



finding of coronal radian ce variation is attributed to a change in total magnetic 



flux in the quiet region. iClose et all (|20G3,) argued that quiet regions are filled 



with short loops connecting magnetic flux fragments in the photosphere, forming 
the so-called magnetic carpet. We speculate that the total magnetic flux of the 
small-scale magnetic loops in the corona changes with the solar activity cycle . 



and hence the heat input to the corona. Indeed, iPauluhn and Solankil (|2003l ) 
studied magnetograms obtained in the quiet regions by SOHO/MDI and found 
that the fraction of network area increased from 10% near the solar activity 
minimum to 13% near the maximum. The long-term variation of magnetic fields 
in the corona needs to be studied in order to understand the coronal heating 
process in quiet regions. Our results should encourage modellers to consider how 
the small-scale magnetic fields in quiet regions are related to the solar activity 
cycle on a global scale. It must be noted that our observations covers only five 
years around the solar minimum. Therefore, the variation amplitude along the 
solar cycle remains to be studied when the solar maximum is reached. 
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